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Abstract: Reflectance spectra measured in Intralipid (IL) close to the 
source are sensitive to wavelength-dependent changes in reduced scattering 
coefficient (/y') and scattering phase function (PF). Experiments and 

simulations were performed using device designs with either single or 
separate optical fibers for delivery and collection of light in varying 
concentrations of IL. Spectral reflectance is not consistently linear with 
varying IL concentration, with PF-dependent effects observed for single 
fiber devices with diameters smaller than ten transport lengths and for 
separate source-detector devices that collected light at less than half of a 
transport length from the source. Similar effects are thought to be seen in 
tissue, limiting the ability to quantitatively compare spectra from different 
devices without compensation. 
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1. Introduction 

Reflectance spectroscopy provides estimates of optical properties within a turbid medium 
such as tissue. This modality has an important role in the field of biomedical optics, as 
reflectance spectra can be used to extract information about aspects of the tissue physiology 
(e.g. vascular oxygen saturation and blood volume fraction) [1,2] and morphology (e.g. cell, 
organelle size and density) [3-6] and can be used to correct raw fluorescence measurements 
for the influences of background scattering and absorption [7-9]. Spectra measured in 
reference media such as the fat emulsion commercialized as Intralipid (IL) are often used to 
calibrate or understand instrument performance [10], yet the spectra recorded by different 
devices can be dramatically affected by device-specific sensitivities to the scattering 
properties. In this study, the influence of detection geometry with respect to IL concentration 
is studied in detail using measurements and Monte Carlo studies, to illustrate and interpret 
how the spectra can change. 

Fiber optic reflectance devices are designed to sample light from a location of interest 
where typically the size and shape is customized for the biological application. Devices that 
collect reflectance remitted at large distances from the source can achieve (relatively) deep 
sampling and therefore are capable of monitoring volume averaged quantities in bulk tissues 
such as blood oxygenation in the brain [11] or breast [12]. Conversely, devices that collect 
reflectance close to the light source sample highly localized volumes, and therefore can 
provide selective sampling of superficial tissue layers [13-19]. Validation of these devices 
with tissue-simulating optical phantoms in vitro is critical to allow characterization of the 
observed reflectance in response to controlled optical properties within the phantom. Phantom 
measurements may also be a valuable method to calibrate the wavelength-specific device 
response for use in analyzing measurements of media with unknown optical properties, such 
as tissue. Perhaps the most commonly adopted scattering phantom is IL [10,20], because it is 
stable across batches, readily purchased and has been extensively investigated to characterize 
the wavelength-dependent optical properties [21-24]. Dilutions of IL can yield wavelength- 
dependent reduced scattering coefficients ( ju[ ) that are similar to values reported in tissue. 

These characteristics have led to wide acceptance of IL as a tissue-simulating scattering 
standard for the investigation of many types of reflectance devices [13,15-19]. However, as 
this study shows, spectra measured in the same IL phantom by different reflectance 
geometries may show substantially different spectral features. 

To understand how scattering properties influence reflectance spectra measured by 
different device geometries differently, it is important to consider the influence of scattering 
on the transport of light being collected during measurements. Light that travels relatively 
large distances in scattering media is termed diffuse, with photon transport occurring in the 
direction of fluence (concentration) gradients within the medium [25,26]. Collected diffuse 
reflectance intensities can be described analytically by the diffusion approximation to the 
radiative transport equation, which allows estimation of both jJ s and the absorption 

coefficient (// fl ) [18,25]. The heuristic guideline for diffuse transport is light that has 
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Fig. 1 . Fiber optic geometries of reflectance devices that measure light close to the source fiber. 
Single fiber measurements obtained using same fiber optic to deliver and collect light with 
diameter given by df. Separate source detector devices have core-to-core separations given by 
rho. Light transport from source to detectors show geometry-specific dependence on the 
probability of scattering events specified by the scattering phase function. Polar plot shows 
phase function for Intralipid at 400 nm. 

travelled 1 to 2 mean free transport lengths from the source [26]; here, the term transport 

length is defined as (fi' s + fi a ) . Analysis of diffuse light is insensitive to the form of the 

scattering phase function (PF), as changes to the scatter frequency (given by ju s ) and 
directionality (given by PF) interchangeably influence diffuse light transport in the form of the 

reduced scattering coefficient ju' s = ju s (1 - g) , where g is the anisotropy given as g = (cos ® s ) 

. As the distance between source and collection location decreases, the sampled transport 
regime shifts from the collection of multiply scattered light to the collection of photons that 
may have undergone few (or single) scattering events [27,28]. Light transport in this semi- to 
non-diffuse regime is not exclusively characterized by //' , but is also influenced by the 
probability of large angle backscattering events, which is defined by the scattering PF [28,29]. 

The influence of PF on light transport collected at various source-detector separations (p) 
is illustrated in Fig. 1 . Light collected at large p is multiply scattered, making it insensitive to 
any individual directional scattering event. Light collected at small p is sensitive to large angle 
scattering events. The frequency of these events is proportional to the magnitude of the 
backscatter-lobe of the PF. In IL phantoms, it is important to consider that both jj' s and PF 
change considerably over the ultraviolet-visible (UV-vis) wavelength range; these changes are 
displayed in Fig. 2 (A) and (B), respectively. While it is well known that fi[ exhibits a power 
law relationship with wavelength [30,31], changes to PF vs. wavelength are less commonly 
understood. In Fig. 2 (B), PFs from selected wavelengths are plotted as the probability density 
distributions of scattering angles; these plots show marked differences between wavelengths 
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Fig. 2. Optical properties of Intralipid optical phantoms. A) shows reduced scattering 
coefficient vs. wavelength for selected IL dilutions. B) scattering phase function of IL at 
selected wavelengths. Right column shows wavelength-dependent metrics of PF: C) 
anisotropy, D) ratio of the first and second moments of the phase function, and the E) 
probability of large angle backscatter events. This study assumes that reduced scattering 
coefficient linearly relates to dilutions of IL, while PF-related parameters are dilution 
independent. 

[24]. Previous studies have developed metrics to describe the average directional scattering 
behavior associated with a PF [28,32]. Specifically, these studies utilized Legendre moments 
of the PF distributions to characterize the average direction of scatter, given by the anisotropy 
( 8 = 8\ )j which is shown in Fig. 2 (C), and the likelihood of large angle backscatter events, 

given by a parameter y = — 



which includes information about the first two moments of 



the PF, g l and g, , respectively and is shown in Fig. 2 (D). For forward-directed scattering 

media, such as IL, y can be interpreted as follows: decreases in values of y are associated with 
an increased likelihood of large angle backscatter events. A simple approach to characterizing 
the differences in the likelihood of backscatter between PFs is to consider the probability of a 
large angle backscatter event, which we define as an scattering angle that would fall within 
the cone of acceptance for the fiber probes included in this study. Consideration of these large 
angle events is especially relevant for the single scattering regime. Figure 2 (E) shows that for 
IL, the likelihood of large angle backscattering events increases substantially across the UV- 
vis wavelength range; correspondingly, y decreases over the same range. Previous studies 
have utilized the moments of the PF to describe collected light intensity near the source 
[28,32-34]. But to the best of the author's knowledge, these findings have not been translated 
to the characterization of scattering-induced spectral features for devices that measure over 
different length scales close to the source, even for a well-characterized scattering standard 
such as IL. 
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This study characterizes scattering features observed in reflectance spectra measured in IL 
optical phantoms. Experimental and simulation data are used to identify and investigate the 
spectral response of reflectance measured in a wide range of IL dilutions as measured by fiber 
optic device geometries with 1) overlapping source and detector with a range of fiber 
diameters, and 2) separate source and detector fibers with a wide range of source-detector 
separations. The underlying mechanisms that link the observed spectral response to the 
wavelength-dependent scattering properties are presented, and the potential errors associated 
with a failure to consider PF-related effects for measurements close to the source are 
discussed. 

2. Methods 

2. 1. Single fiber reflectance experimental setup 

The experimental apparatus for the reflectance device that uses a single fiber optic to deliver 
and collect light during measurements has been described in detail previously [35,36]. The 
device consists of a halogen light source (HL-2000-FHSA, Ocean Optics, Duiven, 
Netherlands), a spectrometer (SD2000, Ocean Optics, Duiven, Netherlands) and a solid core 
single fiber probe that is connected to the light source and the spectrometer via a bifurcated 
fiber. The probe tip is polished at an angle of 15 degrees to remove specular reflections within 
the collected reflectance intensity due to the index of refraction mismatch at the 
probe/medium interface. This study used a calibration procedure to account for variations in 
the output of the lamp, transmission characteristics of the fiber, and sensitivity of the 
spectrometer. Calibration measurements were performed in a calibration phantom containing 
20% IL diluted to a volume particle concentration of 1.5%, which yielded ju' s (800 nm) = 1.2 

mm -1 ( Fjl' (A)), and in a dark container filled with water ( I waler (A) ): 



R° SF {X)=iT{t) 



I phantom (^) ^ water (^) 
jvul 



(1) 



Here, I hantom (A) is the raw signal of a sampled optical phantom, and I s ™ (A) is the absolute 

single fiber reflectance (SFR) signal, expressed in % of incident photons, obtained by Monte 
Carlo (MC) simulations of the IL calibration phantom. Importantly, for these simulations the 
wavelength-specific PFs and //' values reported by [24] were used. The calculation results in 

the calibrated single fiber reflectance intensity, R" SF [X) . 
2.2. Intralipid tissue simulating phantom preparation 

The liquid phantoms were prepared by mixing 0.9% NaCl (Baxter, Utrecht, Netherlands) with 
different amounts of Intralipid 20% (Fresenius Kabi, s-Hertogenbosch, Netherlands); this 
stock solution has a volume particle concentration percentage (IL p ) of 22.7% [22]. Phantoms 
were prepared in volumes of 20 ml with dilutions resulting in IL p values of [0.455, 0.568, 
0.908, 1.135, 1.476, 2.951, 4.540] %. The reduced scattering coefficient of the Intralipid stock 
solution was determined to be //' (800 nm) = 18.6 mm -1 by using the formula given by 

Michels et al. [24]. The dependence of //' on the Intralipid concentration was assumed to be 
linear (for the range of concentrations investigated) and the resulting wavelength-dependent 
ju' s is given for selected phantoms in Fig. 2 (A). The dilutions selected for this study resulted 

in a range of ju' s = [0.4-8] mm -1 that is representative of the values reported in tissue [37]. 

Reflectance measurements were performed by submerging the probe tip a few millimeters 
into the phantom and suspending it multiple centimeters from the container bottom or sides, 
such that the probe would not collect reflections from container surfaces. The measurement 
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was repeated 3 times for each phantom, with these spectra averaged for each phantom. All 
measurements were performed with 5 fiber probes of diameter: [0.2, 0.4, 0.6, 0.8, 1.0] mm. 
The phantom was gently agitated before collecting each spectrum to ensure homogeneity. 
Measurements were reproducible, with variation between repeated measurements of < 3%. 

2.3. Monte Carlo simulations of reflectance measurements 

The Monte Carlo code utilized in this study is based on the steady state light transport in 
multi-layered structures (MCML) approach to stochastically simulate photon propagation 
within a turbid medium [35,38,39]. The code was adapted to simulate the measurement of 
reflectance intensity from a turbid medium with homogeneously distributed optical properties; 
this model of photon propagation has previously been described in detail and will be reviewed 
here briefly. Photons were initialized by selecting a location on the fiber face in contact with 
the medium and were launched into a direction within the fiber cone of acceptance, given as 



9 a = sin 1 , with the fiber numerical aperture (NA) and the index of refraction of the 



medium given by n. Both location and direction were sampled from uniform distributions. 
The fiber probes were modeled in contact with the turbid medium at an air interface, with n of 
the medium set to 1.33 and that of the fiber to 1.45. The fiber NA was specified as 0.22. 
Reflection and refraction at the boundary between fiber face and medium were accounted for 
by using Fresnel's equation and Snell's law. The photons were propagated through the 
medium in discrete steps selected from an exponential distribution that was weighted by the 
scattering coefficient. Each scattering event resulted in a change in propagation direction that 
was sampled from a discrete user-specified PF distribution. Photons that contacted the 
detector fiber within the cone of acceptance were collected. For all simulations, photons that 
propagated far from the fiber face did not contribute to the collected reflectance intensity and 
were terminated at a hemispherical limit from the source fiber at a distance great enough to 
not influence the model outputs for the range of optical properties investigated in this study. 
The collection area for the single fiber simulations overlapped the source fiber location, with 
d f specified in the range [0.2, 0.4, 0.6, 0.8, 1.0] mm. For separate source detector devices, 

source and detector fiber diameters were fixed to 0.2 mm, and photons were collected using 
annular areas at appropriate radial distance such that the core-to-core separation p was 
specified in the range [0.2, 0.4, 0.6, 1.0, 2.0] mm. Simulations returned collected reflectance 
intensity in the absence of absorption, and also returned the weighted average number of 
scattering events experienced by collected photons from each measurement. Reflectance 
intensity returned for separate source detector simulations was calculated by the ratio of the 
total number of photons collected (TPC) by the total number of photons launched (TPL) and 
corrected for differences in collection area (A), as 



The optical properties of the simulated medium were selected to mimic IL phantoms. The 
reduced scattering coefficient was specified by linearly scaling the fi' 3 of the IL 20% stock 

solution, given in [24] for the selected IL dilutions. Simulations of single fiber reflectance 
probes were performed to match the measured optical phantoms, and simulations for the 
separate source detector device were performed for IL dilutions resulting in IL p values of 
[0.454, 0.568, 1.135, 2.270, 3.405, 4.540, 5.675]. Wavelength-dependent PFs were specified 
using empirical formulas also reported in [24], with simulations performed at discrete 
wavelengths between 400 and 800 nm at 100 nm increments. Each simulation initialized 
between 10 and 100 million photons. 





(2) 
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3. Results 



3.1. Reflectance measurement of IL phantoms with a single fiber device 

Figure 3 shows a comparison of reflectance spectra measured in IL optical phantoms by 
reflectance devices that use a single fiber optic to deliver and collect light. Panels (A), (B) and 
(C) contain spectra measured in a dilution series of IL phantoms, with IL p in the range [0.454- 
4.54]%, with each panel showing spectra obtained by a different fiber diameter from the series 
[0.2, 0.4, 1.0] mm. Initial inspection of the data show expected trends of increased R" SF (X) in 
response to increases in IL p , which is observed by comparison of spectra within each panel. 
R° SF {X) also increases with increases in d f , which is observed by comparison of spectra 
between different panels. A closer inspection of the spectra shows that the wavelength- 
dependent spectral shape of each curve changes with respect to changes in IL p and d f . Figure 

3 (D) and (E) show normalized R° SF (/l) spectra measured by 3 different fiber diameters in 

selected IL p phantoms of 4.54% and 0.454%, respectively. Spectra are normalized to the 
intensity at 400 nm to highlight the substantial differences in spectral shapes returned by 
measurements with different fiber diameters. First we consider the spectra measured by d f = 
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Fig. 3. Analysis of spectral scattering features observed in experimentally measured reflectance 
collected by single fiber reflectance device. Left hand column presents calibrated reflectance 
spectra as measured with single fiber diameters of A) 0.2 (x marks), B) 0.4 (circle marks), and 
C) 1.0 mm (square marks). Each panel shows spectra measured from dilutions of IL resulting 
in volume particle concentrations of [0.455, 0.568, 0.908, 1.135,1 .476, 2.95 1, 4.540]% that are 
indicated by colors [red, blue, green, yellow, magenta, cyan, and black], respectively. Middle 
column shows selected spectra measured by each fiber in D) 4.54% and E) 0.454%, with each 
spectra normalized at 400 nm for comparative purposes. Right hand column shows the average 
number of scattering events experienced by photons at discrete wavelengths for spectra in F) 
4.54% and G) 0.454% phantoms. 
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1 mm in the ILp = 4.54% phantom, where R'^ F {X) monotonically decreases vs. wavelength. 
Reflectance spectra measured by smaller fibers in the same phantom show a transition to a 
more complicated biphasic behavior. Specifically for d f = 0.2 mm, R" SF {A) monotonically 

decreases over 400-500 nm, but exhibits a plateauing effect over 500-800 nm. For the lower 
scattering phantom, with IL p = 0.454%, reflectance spectra exhibit the biphasic behavior for 
d f = [0.2, 0.4, 1.0] mm, with R° F (A) decreasing over 400-500 nm followed by R^W 
increasing from 500 to 800 nm. The data in panels (D) and (E) show that the wavelength- 
dependent changes in R° F measured in the same phantom are substantially different between 
different fiber diameters. For example, reflectance spectra measured in the 4.54% phantom 
show a decrease in R° F from 400 to 800 nm of 25% for d f = 0.2 mm, and 45% for d f = 1.0 

mm. Moreover, the R" SF data show substantially different spectral responses for measurements 
of phantoms with different IL dilutions with the same fiber diameter. For example, reflectance 
spectra measured with a fiber with d f = 0.4 mm in the 4.54 and 0.454% IL p phantoms show 

of decrease of 40% and an increase of 5%, respectively. 

Understanding the mechanisms underlying these spectral differences requires 
consideration of two contributing factors: 1) how scattering influences light propagation 
during measurement and, 2) how these influences change over the wavelength range sampled 
in the IL phantoms. It is understood that reflectance devices sample photons that have 
experienced a range of scattering events [27]. For reflectance measured near the source, 
decreases in /u' s or d f are likely to decrease the number of scattering events experienced by 
collected photons [40], which in turn increases the sensitivity of the collected reflectance 
intensity to PF-effects. In IL the scattering properties that influence light transport ( )J s and PF) 

change substantially across the UV-vis wavelength range. Figure 2 (A) shows n' s (X) over 

400-800 nm. The PF also changes over the same wavelength range, with y{X) showing that 

the likelihood of backscatter events increases vs. wavelength. To characterize the link 
between these wavelength-dependent optical properties and the observed reflectance, we 
considered the scattering histories of photons collected across the sampled wavelength range. 
For this purpose, a Monte Carlo model was used to return the average number of scattering 
events for photons collected at discrete wavelengths along each spectrum reported Fig. 3 (D) 
and (E); the corresponding scattering event data are shown in panels (F) and (G), respectively. 
These scattering data are plotted on a log-scale and include annotations to highlight two 
important thresholds: 1) the single scattering baseline, and 2) an approximate limit to the 
diffuse regime of 10 scattering events. We first consider the reflectance spectrum that 
exhibited single -phase behavior, (d f = 1 mm measurement of the 4.54% IL p phantom). Light 

collected across this spectrum was, on average, multiply scattered (with >10). The resulting 
diffuse signal is expected to be insensitive to PF-effects, a claim that is supported by 
agreement in wavelength-dependent trends of R'^ F (A) and // s '(/t) over the sampled 

wavelength range. Inspection of scattering data for reflectance curves that exhibit the convex 
biphasic behavior shows that the collected photons often experienced few (<10) scattering 
events. These measurements contained contributions from photons sampled in the semi- 
diffuse regime and are expected to be sensitive to the wavelength-dependent changes in PF. 
This claim is supported by considering that the observed increase in R^ F (A) vs. wavelength 
in the 500-800 nm region for these spectra is correlated with the wavelength-dependent 
changes in backscatter probability described by y(X) and PF( 0 a ) as observed in Fig. 2. These 

data indicate that the spectral features observed in R" SF (A) spectra measured in IL phantoms 
are attributable to the sensitivity of the reflectance device to the combined wavelength 
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dependent changes in fj! s (X) and PF in IL. While the combined influence of n\ and PF on 
R" SF at a single wavelength has been modeled in detail previously [34,41]; these data present a 
qualitative representation of how wavelength-dependent features in reflectance spectra are 
affected by the wavelength dependent scattering properties and device geometry. 

The spectral features observed in single fiber reflectance spectra measured in IL phantoms 
highlight the changing influence of PF across the wavelengths of a single spectrum; these 
changes can have profound effects on the relationship between R" SF and ju' s . The left column 

of Fig. 4 shows R^ F vs. fj! s measured in all IL phantoms with panels (A), (B), and (C), 
showing measurements by d , = [0.2, 0.4, 1.0] mm, respectively. Here, colors are used to 
indicate wavelength. These data show that R° SF cannot be used to uniquely identify fi[ , 
because PF-effects influence the proportionality between Rg F (A) and fi' s differently at each 
wavelength. This statement is evidenced by the data in Figs. 3 (D) and 3(E), where R^W 
does not explicitly follow n' s (X) , but is also influenced by y{A) . Comparison of 
measurements from different fiber diameters showed an underlying relationship between R° F 
and the dimensionless reduced scattering coefficient, given as ju' s d f . Figure 4 panel (D) 
contains measurements from all investigated fiber diameters, d f [0.2, 0.4, 0.6, 0.8, 1.0] mm, 
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Fig. 4. Relationship between experimentally measured single fiber reflectance intensity and 
scattering properties in IL phantoms. Left hand column shows reflectance intensities vs. 
reduced scattering coefficient for fiber diameters of A) 0.2, B) 0.4, and C) 1.0 mm. Panel D) 
shows reflectance intensities vs. dimensionless reduced scattering for all measured fibers. In 
this plot color indicates different wavelengths and can be used to interpret the phase function 
specific proportionalities between reflectance intensity and reduced scattering coefficient, with 
colors transitioning from blue (at 400 nm) to red (at 800 nm). 
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in investigated phantoms with IL p in the range of [0.454-4.54]%. Inspection of the R" SF data 

show that for increasing fj,' s d f values, less stratification is observed between different 

wavelengths, meaning that the collected intensities are less influenced by PF. This is 
consistent with previous work that identified a threshold limit of PF-independent behavior for 

values of ju' s d f >10 [35]. This plot shows that reflectance intensity collected by a single fiber 
optic is interchangeably influenced by jf s and d f and highlights the influence of the 
wavelength-dependent PF on R° SF . 
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Fig. 5. Analysis of spectral scattering features observed in simulated reflectance spectra 
collected by devices with separate source and detector. Left hand column shows dimensionless 
reflectance spectra returned using selected source detector separations of A) 0.2 (x marks), B) 
0.4 (circle marks), and C) 1.0 mm (square marks). Each panel shows spectra simulated from 
dilutions of IL resulting in volume particle concentrations of [0.455, 0.568, 0.908, 1.135, 1.476, 
2.951, 4.540]% that are indicated by colors [blue, red, green, yellow, magenta, cyan, and 
black], respectively. Middle column shows spectra returned by 3 selected source-detector 
separations in D) 4.54% and E) 0.454% phantoms, with the spectra normalized at 400 nm for 
comparative purposes. Right hand column shows the average number of scattering events 
experienced by photons at discrete wavelengths for spectra in F) 4.54% and G) 0.454% 
phantoms. 

3.2. Reflectance measurement of IL phantoms with separate source-detector fibers 

Figure 5 shows simulated reflectance measurements in IL optical phantoms performed with 
separate source and detector fibers. Here, panels (A), (B), and (C) contain spectra measured in 
a dilution series of IL phantoms, in the range [0.454-5.675]%, with each panel showing data 
obtained by a source-detector separation from the series p = [0.2, 0.4, 2.0] mm. Initial 
inspection of these data show a complicated dependence of R^ D (A) on the phantom IL p . For 
adjacent source-detector fibers (p = 0.2 mm), an increase in IL p results in an increased 
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reflectance intensity across the entire spectrum. This is not the case for source detectors with 
large separation (p = 2.0 mm), as increases in IL p from 0.4 to 1.0% result in increased R^ D , 

but increases in IL p from 2 to 5% result in decreased Rg D ; this nonlinear dependence of R^ D 
on n' s is described by diffusion theory and has been characterized in detail previously [18]. In 
order to understand the influence of scattering properties on the spectral response, we first 
compare spectra measured in the same phantom with different p. Figure 5 panels (D) and (E) 
show reflectance spectra measured at 3 different source detector separations in phantoms with 
IL p of 4.54% and 0.454%. The panels (F) and (G) show the average number of scattering 
events corresponding to discrete wavelengths in the spectra plotted on panels (D) and (E), 
respectively. The normalized spectra from the 4.54% IL p phantom show decreases in R^ D (X) 
vs. wavelength for the small fiber separations (p= [0.2, 0.4] mm), but an increase of 

vs. wavelength for the large p (2.0 mm). The latter spectrum can be described by increases in 
ju' s vs. wavelength increasingly attenuating the collected reflectance intensity; therefore 

R^,(A) follows the inverse trend of fJ s {X) vs. wavelength over the sampled wavelength 
range. Inspection of the scattering data for these spectra in panel (F) show that these 
measurements sample diffuse light, and are independent of PF. Figure 5 panel (E) shows 
normalized reflectance spectra measured in the 0.454% IL p phantom at different p. These data 
show a transition from Rg D {X) monotonically decreasing vs. wavelength (for p= 2.0 mm) to 

biphasic behavior (for p= 0.2). Inspection of the number of scattering events experienced by 
collected photons in panel (G), shows a transition from sampling multiply scattered light in 
the diffuse regime (for p = 2.0 mm) to the sampling of fewer scattering events in the semi- 
diffuse regime (for p = 0.2 mm). These observations are consistent with the concept that 
measurements close to the source are more likely to be sensitive to the frequency of large 
angle backscatter events, and therefore the collected spectrum is influenced by PF-effects. 
Figure 6 shows analysis of the relationship between R^ D and fi' s for separate source 

detector devices. Panels (A), (B), and (C) show R^ D vs. /j' s for measurements of phantoms 
with a wide range of IL p , with each panel showing measurements from a device with p = [0.2, 
0.4, 2.0] mm. R^ D measurements close to the source (p = 0.2 and 0.4 mm) show a sensitivity 
to PF effects for small /J values; this is characterized by the non-unique relationship between 
R^ D and //' in the low scattering region. For p = 0.2 and 0.4 mm devices, as ft' increases the 
R^ D shows less stratification due to PF-effects and the data transition to PF-insensitivity for 
large fi' s values; this transition follows the shift from sampling predominantly in the semi- 
diffuse regime to the diffuse regime as scattering increases. The data for the large source 
detector separation of p = 2.0 mm shows a distinct nonlinear dependence on fj! s , with the 

Rg D vs. n\ showing a monotonically increasing region that transitions through a peak into a 
monotonically decreasing region; this transition occurs for measurements sampling diffuse 
light that are independent of PF effects. Inspection of the reflectance data measured at 
different p shows an underlying relationship that is identified by conversion of the x-axis to 
the dimensionless reduced scattering coefficient, here given as /u' s p , and conversion of the y- 

axis to Rs D p 2 ■ Figure 6 panel (D) contains measurements from all investigated source- 
detector separations/? = [0.2, 0.4, 0.6, 1.0, 2.0] mm in all simulated IL phantoms IL p = [0.454- 
5.675] %. These data show that Rs D p 2 is influenced similarly by fi' s p for different p. 

Inspection of these data showed that for ju' s p >0.5 , R* D showed less than 10% variation 
about the mean value, and are therefore considered independent of PF influences at different 
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Fig. 6. Relationship between reflectance intensity returned by simulations of separate source- 
detector geometries and scattering properties in IL phantoms. Left hand column shows 
reflectance intensities vs. reduced scattering coefficient for source-detector separations of A) 
0.2, B) 0.4, and C) 2.0 mm Panel D) shows reflectance intensities vs. dimensionless reduced 
scattering for all measured separations. In this plot color indicates different wavelengths with 
[400, 500, 600, 700, 800] nm indicated by [blue, green yellow, black, red], respectively. 
Vertical dashed lines indicate the empirically identified threshold for observed PF-sensitivity of 
the reflectance at small dimensionless scattering values. 

wavelengths. This threshold was used to identify the smallest /u' s that could be sampled 
independent of PF for each p. For separate source detectors of p = [0.2, 0.4, 2.0] mm, the 
lower bound of /j' s was estimated as [2.5, 1.25, 0.25] mm -1 , respectively. These thresholds are 

included as vertical dashed lines on plots in Fig. 6 to allow visualization of scattering regions 
influenced by PF and independent of PF. 

4. Discussion 

IL optical phantoms have substantial wavelength-dependent changes in p! s and PF in the UV- 

vis region. This study characterizes how these scattering properties influence spectral features 
of reflectance measured close to the source. Experiments and simulations were used to 
investigate reflectance spectra returned in a set of IL phantoms that contained a wide range of 
fi[ values. Reflectance spectra were measured or simulated for device geometries that 
delivered and collected with either a single fiber or separate fibers. Data show that changes to 
the observed spectral remission characteristics occur in response to slight changes in either 
measurement geometry or IL dilution. These changes are related to the light transport regime 
sampled, with spectral features showing a (non-linear) dependence on fi[ for measurements 

that collect diffusely scattered light, and on top of that a dependence on PF for measurements 
that collect semi-diffuse light. 
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Tissue simulating optical phantoms are commonly used for two purposes: 1) 
characterization of observed reflectance response to controlled optical properties within the 
phantom, and 2) calibration of the wavelength-specific reflectance device response to allow 
system independent interpretation of measurements made in samples with uncontrolled optical 
properties. The results presented in this study have implications for both of these uses. Device 
characterization studies may seek to establish or validate a relationship between observed 
reflectance and ju' s [16-18]. However, the data presented in this study show that scattering PF 
can profoundly influence observed relationships between reflectance intensity and fj! s . The 
data in Figs. 4 and 6 show that for reflectance measurements that sample dimensionless 
reduced scattering coefficients below empirically identified threshold values (n' s d f =10 for 

single fiber devices and /j! s p = 0.5 for separate source detector devices) the collected 

reflectance intensities are sensitive to PF-effects. Reflectance measurements performed below 
these thresholds collect a substantial component of light that has undergone few, or single, 
scattering events, making the probability of collecting light in that regime (at least somewhat) 
dependent on the backscattering features of the scattering PF. These data indicate that a 
proportionality established between reflectance intensity and ju' s in IL at a single wavelength 

will not hold at other wavelengths. Moreover, proportionalities are not directly extensible to 
another scattering medium that may contain an unknown wavelength-dependent PF, such as 
tissue. 

IL optical phantoms are also commonly used for system calibration purposes [10]. While 
this approach can be valid, care must be taken to consider influences of the PF. As shown in 
Figs. 3 and 5, spectral remission features can depend strongly on the IL dilution of the 
phantom. It is possible to account for these features by characterizing the response of the 
investigated device to the moments of the PF. This can be investigated either by experimental 
measurements in another scattering medium with a controlled wavelength-dependent PF, such 
as a phantom constructed with polystyrene beads [42], or simulation, where variation of PF 
can be completely controlled. The experimental data for single fiber reflectance spectra in this 
study were calibrated with measurements in IL, by using a Monte Carlo analysis to convert 
the measured reflectance spectra to an absolute scale (in units of % of incident photons) [35], 
a conversion that accounts for the wavelength-dependent PF influence. Such a calibration 
becomes independent of IL p This approach to spectral analysis requires sufficient information 
to estimate PF effects; such information could be provided by spectral information measured 
at an additional length scale such as with multiple fiber diameters [34,35,41] or at multiple 
distances from the source [14,33], or spectral information measured with an angular resolution 
[43,44]. The use of IL as a calibration phantom without correction for PF-effects may 
introduce a substantial bias into quantitative estimates of scattering parameters, especially if 
applied to measurements in a sample of unknown optical properties, such as tissue. 

This study identified an interesting difference in the relationship between reflectance 
intensity and fi' s for device geometries that utilize a single fiber as a source and detector 
compared with a device that utilizes separate source and detector fibers. Reflectance intensity 
collected by single fiber devices shows a monotonically increasing response to increases in 
ju' s . , as displayed in Fig. 4. In contrast, reflectance intensity collected by separate source 

detectors shows a nonlinear response that increases or decreases based on fi' s p , as shown in 

Fig. 6. These limiting behaviors can be intuitively explained. The single fiber geometry allows 
collection of light that has scattered an infinitesimally small distance into the medium prior to 
collection. Therefore, in the limit of highly scattering media, the single fiber reflectance 
device is no less likely to collect light in response to increased fj! s . In contrast, the collection 
geometry with separate source and detector fibers has a geometrical minimum path that light 
must travel within the medium prior to collection. Therefore, in the limit of highly scattering 
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media, the collected signal will be attenuated in response to increases in ju' s . These 
differences in the required light path travelled by light collected by either device leads to 
different threshold sensitivities to PF, which empirically are fi' s d f = 10 and jJ s p = 0.5 for 

each device-type, respectively. This means that for single fiber measurements, the transport 
length, defined by the inverse of the reduced scattering coefficient, must be an order of 
magnitude greater than the collection length to be independent of PF effects. A device with 
separate source detectors is (in general) less sensitive to PF, with collection at locations up to 
half of the transport length showing insensitivity to PF. This threshold is expected to be 
dependent on the source and detector fiber diameters, which were fixed at 0.2 mm in the 
current study; the exact value of the threshold may vary for different fibers sizes. It should be 
noted that in IL phantoms constructed to mimic the range of fj! s expected in tissue, in the 

range of [0.5-2.0] mm -1 [37], devices that collect with p as large as 1 mm can be influenced 
by PF effects. Moreover, inspection of the data in Figs. 5 and 6 show that measurements made 
in such phantoms with adjacent source and detector fibers (p = 0.2 mm) may contain a 
spectral signature that shows little to no wavelength-dependence on jJ s . 

While this study specifically investigates the spectral response of specific reflectance 
geometries, the findings are important for other measurement devices as well. This study 
characterized the response over a set range of geometrical perturbations to single fiber size 
and source-detector distance. This study does not perform a comprehensive evaluation of all 
possible geometrical configurations of fiber probes (e.g. all possible source and detector fiber 
sizes, inter- fiber distances, or numerical apertures). The results of this study are applicable to 
reflectance spectra measured by any probe configuration: the sampled transport regime will 
determine sensitivity to wavelength-dependent changes in PF. The results are also important 
for devices that sample reflectance at a single wavelength. In this case, care must be taken to 
correctly interpret the collected signal to avoid mischaracterizing the observed reflectance 
intensity as influenced exclusively by (j! s . For example, the use of raw reflectance intensity at 
a single wavelength to correct fluorescence for the influence of optical properties will be 
problematic considering that fi' s and PF will influence reflectance and fluorescence intensities 

differently [9]. Ongoing work will characterize the PF-specific influences of fluorescence 
correction algorithms. 

5 Conclusions 

This study investigates the influence of the wavelength-dependent scattering properties /j! s 

and PF on spectral features of reflectance collected near the source location in IL optical 
phantoms. Experiments and simulations are used to characterize features in the 400-800 nm 
wavelength range from IL optical phantoms with IL p in the range 0.454-5.675%, which relates 
to a wide range of ju' s = [0.3-8] mm -1 . Data were measured or simulated for reflectance 
geometries that 1) utilized a single fiber optic for both source and detector with the d, = 
[0.2-1.0] mm, and 2) separate source and detector fibers with core to core separations in the 
range p = [0.2-2.0] mm. Results for both collection configurations showed a dependence of 
spectral scattering features on the transport regime sampled during measurement. PF- 
dependent effects were observed in reflectance spectra measured below a dimensionless 
threshold of fi' s d f =10 (i.e. ten transport lengths) for single fiber devices and fJ s p = 0.5 (i.e. 

half a transport length) for separate source detector devices. The results show the need to 
consider the influence of PF when using measurements of IL to either characterize or calibrate 
a reflectance device that collects light close to the source. 

Acknowledgments 

This work has been funded by NIH grants POl CA084203 and POl CA080139. 



#163835 - $15.00 USD Received 2 Mar 2012; revised 12 Apr 2012; accepted 20 Apr 2012; published 24 Apr 2012 
(C) 2012 OSA 1 May 2012 / Vol. 3, No. 5 / BIOMEDICAL OPTICS EXPRESS 1 100 



